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Multiple parton scattering and induced parton energy loss are studied in deeply in- 
elastic scattering (DIS) off nuclei. The effect of multiple scattering of a highly off-shell 
quark and the induced parton energy loss is expressed in terms of the modification 
to the quark fragmentation functions. We derive such modified quark fragmentation 
functions and their QCD evolution equations in DIS using the generalized factor- 
ization of higher twist parton distributions. We consider double-hard and hard-soft 
parton scattering as well as their interferences in the same framework. The final 
result, which depends on both the diagonal and off-diagonal twist-four parton distri- 
butions in nuclei, demonstrates clearly the Landau-Pomeranchuk-Migdal interference 
features and predicts a unique nuclear modification of the quark fragmentation func- 
tions. 



I. INTRODUCTION 

Hard processes in high-energy strong interactions are always localized in space-time because of 
the large momentum-energy transfer. The asymptotic behavior of QCD allows one to compute 
these cross sections perturbatively. Together with the factorization theorem and the experimental 
information of parton distributions and fragmentation functions, hard processes in hadronic collisions 
have been well understood |Q. One can then in turn use them as probes of nuclear matter as well 
as hot quark-gluon plasma which is expected to be formed in high-energy nuclear collisions. In 
particular, the propagation of an energetic parton and its induced energy loss has been proposed as 
a probe of the properties of dense matter formed in high-energy nuclear collisions ■ Based on a 
model of multiple scattering and induced radiation in QCD proposed by Gyulassy and Wang (GW) 
H , recent theoretical studies show that a fast parton will lose a significant amount of energy via 
induced radiation when it propagates through a hot partonic matter. The most interesting feature of 
the result is the quadratical distance dependence of the total energy loss because of the non-Abelian 
nature of QCD radiation and the Landau-Pomercanchuk-Migdal (LPM) jl(J interference. Such a 
quadratical distance dependence is also a consequence of the GW static color-screened potential 
model for multiple scattering where colors are screened but not confined. In this paper we will study 
parton multiple scattering inside a nucleus where colors are confined to the size of a nucleon. In this 
case parton propagation will certainly be different from that in a partonic matter and one should 
expect that the parton energy loss to be related to the nucleon size or confinement scale. Parton 
energy loss in eA DIS has been studied before within various models of intranuclear scattering and 
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the modification of hadronization |llj |15[| . In this paper we will study the problem in the framework 
of multiple parton scattering in perturbative QCD (pQCD). 

Unlike the situation in QED, the energy loss of a parton cannot be directly measured because 
partons are not the final experimentally observed particles. The total energy of a jet as traditionally 
defined by a cluster of hadrons in the phase space will not change much due to medium induced 
radiation because a jet so defined contains particles both from the leading parton and from the 
radiated gluons. This is particularly the case if multiple scattering and induced radiation do not 
dramatically change the energy profile of the jet in phase-space. It is also virtually impossible 
to determine the jet energy event by event because of the large background and its fluctuation 
in heavy-ion collisions. One then has to resort to particle distributions within a jet and study 
the effect of parton energy loss by measuring the modification of the particle distribution due to 
multiple scattering and induced radiation ||. One such distribution is the fragmentation function 
of the produced parton, D a ^h(z, /x 2 ), where z is the fractional energy of the parton a carried by the 
produced particles h. Unlike the situation that has been considered in most of the recent theoretical 
studies of parton energy loss of an on-shell parton during its propagation through medium, 

partons produced in hard processes are normally far off-shell characterized by the momentum scale 
y 2 . Final-state radiation of these off-shell partons in free space, such as in e + e~ annihilation, leads 
to the ^(-dependence of the fragmentation functions as given by Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) @ QCD evolution equations, 

/ — [jq^qg{z)D q ^ h {zh/z,H ) + ~f q ^g q {z)D g ^ h (z h /z, fl )J , (1) 
Jz h z 



ain M 2 2tt 
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lg^qq( Z ) D q^h(Zh/z, H % ) + Jg^gg (z)D g ^ h (Z h /Z, [J?) 

9=1 



(2) 



where j a —>bc(y) are the splitting functions of the corresponding radiative processes fll7| . When a 
parton is produced in a medium, it will suffer multiple scattering and induced radiation that will 
then lead to medium modification of the DGLAP evolution of the parton fragmentation functions. 
In this paper we will derive such modified evolution equations for parton fragmentation functions in 
the simplest case of deeply inelastic eA scattering (DIS). Multiple scattering and induced radiation 
suffered by the leading quark in this case will give rise to an additional term in the DGLAP evolution 
equations. As a consequence, the modified fragmentation functions become softer. This can be 
directly translated into the energy loss of the leading quark. 

The study here is very similar to that by Luo, Qiu and Sterman (LQS) |l8[ on nuclear depen- 
dence of jet cross section in deeply inelastic scattering. The difference is that they considered large 
transverse momentum jet production whereas we will concentrate on soft gluon emission which is 
responsible for the DGLAP evolution equations of fragmentation functions. Depending on the frac- 
tional momentum carried by the second parton in the case of double parton scattering, one can 
categorize the processes into soft-hard, hard-hard and their corresponding interferences. In their 
study of large transverse momentum jet production, LQS only considered soft-hard and hard-hard 
processes, because the interference terms become negligible as we will show. Since we have to 
consider soft gluon emission, the interference terms are very important and they will cancel contri- 
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butions from hard-hard and hard-soft processes in the limit of collinear emission (or zero transverse 
momentum). In this paper, we will include all these processes and treat them in the same manner. 
Utilizing the generalized factorization of higher-twist parton distributions, we find that each process 
in the double scattering probes different twist-four parton correlation of the nuclear medium. As a 
result of the sum of all contributions, the additional term in the modified DGLAP evolution equation 
is proportional to the combined twist-four matrix element, 

J dy'dy^dy^e(-y^)9(y- - y^)(l - e **M> + (»r-ir))(i - e ~ ix ^) 

^ B+XL ) P +y- (A\^ q (0)~f+F a +(y-)F+(y-)Tjj q (y-)\A), (3) 

where 

x L = ^/2p+q-z{l- z), (4) 

with [i the typical transverse momentum or factorization scale and z the fractional momentum of 
the emitted gluon. In this paper, we denote the four-momentum of the virtual photon and the target 
nuclcon in DIS as 

q=[-Q 2 /2q-,q-,6 1 _] 1 

p=[p+,0,6 ± l (5) 

respectively. The Bjorken variable is then = Q 2 /2p + q~ . 

In the above matrix element, one can identify 1/xlP + = 2q~ z(\ — z)/ [i 2 as the formation time 
of the emitted gluons. For large formation time as compared to the nuclear size, the above matrix 
element vanishes, demonstrating a typical LPM interference effect. This is because the emitted gluon 
(with long formation time) and the leading quark are still a coherent system when they propagate 
through the nucleus. Additional scattering will not induce more gluon radiation, thus limiting the 
energy loss of the leading quark. 

The reduction of phase space available for gluon radiation due to the LPM interference effect is 
critical for applying the LQS formalism to the problem in this paper. In the original LQS approach 
[pi, the generalized factorization for processes with a large final transverse momentum t\ ~ Q 2 
allows one to consider the leading contribution in 1/Q 2 which is enhanced by the nuclear size 
Ra ~ A 1 / 3 . For large Q 2 and A, higher- twist contribution from double parton rescattering that is 
proportional to ol s RaIQ 2 will then be the leading nuclear correction. One can neglect contributions 
from more than two parton rescattering. In deriving the modified fragmentation functions, we 
however have to take the leading logarithmic approximation in the limit dp <C Q 2 , where £t is the 
transverse momentum of the radiated gluon. Since the LPM interference suppresses gluon radiation 
whose formation time (r/ ~ Q 2 /^tP + ) is larger than the nuclear size MRa/p + in our chosen frame, 
t\ should then have a minimum value of t\ ~ Q 2 /MRa ~ Q 2 /A 1 ^ 3 . Here M is the nucleon mass. 
Therefore, the logarithmic approximation is still valid for large nuclei (MRa 3> 1). In the meantime, 
the leading higher-twist contribution proportional to a s RAji\ ~ cc s Ra/Q 2 can still be small for 
large Q 2 so that one can neglect processes with more than two parton rescattering. The parameter 
for twist expansion of the fragmentation processes inside a nucleus is thus a s A 2 / 3 /Q 2 as compared 



3 



to ctsA 1 / 3 /Q 2 for processes with large final transverse momentum as studied by LQS ]l8| . This 
is why the nuclear modification to the fragmentation function as derived in this paper depends 
quadratically on the nuclear size Ra- 

Since our approach allows us to relate the energy loss of a fast parton to the parton correlations 
in the medium, the final results will be sensitive to the properties of the medium through which a 
produced parton propagates. In particular, the medium-dependence of the above matrix elements 
in a deconfincd quark-gluon plasma will be very different from that in an ordinary nuclear medium. 
Therefore, the parton energy loss and its dependence on the medium size in hadronic matter are 
expected to differ from that in a quark-gluon plasma. This simply reflects the different parton 
correlations in different types of media. 

The results of the present study were already reported in Ref. |^9| . We will provide a detailed 
derivation and discussion in this paper. The rest of this paper is organized as follows: In the next 
section, we give a brief overview of the framework of our study including generalized factorization 
of twist-four parton distributions in hard processes. In Section III we will discuss in detail the 
calculation of different contributions to soft gluon emission involving double parton scattering. To 
simplify the calculation, we will use the techniques of hclicity amplitude with soft gluon approxi- 
mation (z << 1). In Section IV, we will consider virtual corrections from unitary constraints. We 
will then define, in Section V, the effective parton fragmentation functions in deeply inelastic eA 
collisions and derive the modified DGLAP evolution equations. We will also calculate the average 
energy loss suffered by the leading quark. Finally, in Section VI, we will discuss our results and their 
implications in other hard processes. 



II. GENERAL FORMALISM 



Consider the following semi-inclusive process in the deeply inelastic lepton-nucleus scattering, 

e(Li) + A(p) — » e(L 2 ) + h{l h ) + X , (6) 

where L\ and Li are the four-momenta of the incoming and the outgoing leptons respectively, p 
is the momentum per nucleon for the nucleus with the atomic number A, and £h is the observed 
hadron momentum. The momentum of the virtual photon (7*) is q = L2 — L\. 

The differential cross section of semi-inclusive processes in DIS with an observed final state hadron 
lh can be expressed as 

where s = (p + L\) 2 is the total invariant mass of the lepton- nucleon system. The leptonic tensor 
L M „ is given by 

= i Tr(7 • Li7 p 7 • L 2 J U ) ■ (8) 
The semi-inclusive hadronic tensor Ei h dW^ u / 'cPi^ is defined as, 
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Ei h -^j- = - J2( A \ J/1 (°)l^ M*' h\r(0)\A)2nS 4 (q+p- p x 4) , (9) 

h x 

where J^x runs over au possible intermediate states and is the hadronic electromagnetic (EM) 
current, J M = e^^Vv 




FIG. 1. Diagram representing the factorized form of the semi-inclusive hadronic tensor Wpu. 

In the parton model with collinear factorization approximation, one can in general factorize the 
semi-inclusive cross section into a parton fragmentation function and the partonic cross section. By 
factorizing the parton fields out of the EM current as 

j"(o) ^£j£ Q (o)C(o) = £^(o)j£t (o) , (io) 

c,a c,a 

one can express the semi-inclusive hadronic tensor as 

= £ Z^Jz,,4)E^I J ^( )l^)^l J ^(°)l^) 2 ^c) , (11) 



dzh 

as illustrated by Fig. [l], where a, ft are indices (spinor for quarks and Lorentz for gluons) of the 
partonic field of species c with momentum £ c , and d"z^ h (zh, &h), related to the fragmentation function, 
is defined as 

%t h (z h ,£ h ) = £/ ^dMO\k(0)\h,S){h,S\^(y)\0}e-^yS(z h - ^). (12) 

After such factorization, one only needs to calculate the cross sections of partonic hard processes. 
The fragmentation functions defined in terms of parton matrix elements in the above equation are 
nonperturbative and can only be obtained via experimental measurements. In addition, one can 
also factorize out the parton distributions inside nuclei (or nucleons). Such factorization has been 
proven to all orders and in the leading twist . 

In the infinite- momentum frame as defined by Eq. (^), the dominant momentum component of 
the leading parton in a single jet event is the minus component. Therefore, we define the momentum 
fraction carried by a produced hadron as Zh — lZ/£~. 
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A. Leading Twist Contributions 
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FIG. 2. Lowest order process that contributes to H^O),. 



The leading- twist contribution to DIS to the lowest order (0(a°)) comes from a single hard 7* + q 
scattering as illustrated in Fig. ||. The semi-inclusive hadronic tensor dW^/dzh can be factorized 
as M 



dw m r 



(13) 



where f q {x) is the quark distribution as defined by 



^I e ^-i(A|^(0) 7 + ^(y-)|A). 

Z7T Z 



Here quark fragmentation function D q ^h{zh) is defined as 



(14) 



D q ^ h (z h ) = -fd q ^ h (z h ,£ h ) = -pTr[-y-d q ^ h (z h ,e h )], 
d a q i h {z h J h )=Y, J ^d i ye-^yS(z h -^)(0\^ q (0)\h,S)(h,S\^(y)\0) 

,l!,+ e-^y^^(0\^ q (0)\h,S)(h,S\^(y + )\0) . 



z, 2 / 2tt 



The hard part of 7* + q partonic scattering is 



^ = < 2 Tr(7 ' ' (g + (2?r),5[(<z + ^ 



l i r T 
2 



(5(a; — xb) 



(15) 



(16) 



(17) 



where momentum conservation gives £ q = xp+q and by definition £ h = Z] l l q = z^q ■ The transverse 
and longitudinal tensors are defined as 



,L 



1 r p- q ir p- 9 , 



(18) 
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Note that a sum and an average over color indices are implied in the quark distribution and fragmen- 
tation function, respectively. In Eq. (|l5|), we used the collinear approximation in the unpolarized 
fragmentation functions, 



(19) 



while other higher-twist terms are neglected. We will use such a collinear approximation throughout 
this paper even in the case of double scattering. We will also neglect other higher-twist contributions 
to the fragmentation function that are independent of the nuclear size. 




FIG. 3. The hard partonic part of next-to-leading order process that contributes to H$ . 



At the next-to- leading order 0(a s ), the dominant (in leading log approximation) real radiative 
contribution to the fragmentation process in an axial gauge (A~ = 0) comes from the hnal state 
radiation. Since we are only interested in the collinear behavior of the radiative corrections in 
order to study the evolution equation of the fragmentation functions, we will keep only the leading 
(divergent) contribution when the gluon's transverse momentum vanishes(with respect to q + xp), 
i.e., £t — * 0. Using the same collinear approximation, one finds the leading contribution from the 
radiative correction (for £j, up to a factorization scale /j 2 ) to the quark fragmentation process, 



dW, 



S(l)q 



I dxf A (x) f 1 -D q ^ h (z h /z)H^(x,p,q,z). 

„ J J Zh. Z 



dz h 



(20) 



where Zh = ^h/Q^ an( l z = jq~ are the fractional momentum carried by hadron and the final 
quark, respectively. The factorized form in Eq. ( ^p|) is quite general for all processes involving gluon 
radiation in the final state. It is then enough to calculate the partonic hard part H^) q (x,p, q, z) as 
illustrated in Fig. ||, 



H$«(x,p, q ,z)=Hj?J(x, P , q )£ j£%C F \ 



1 + z 2 



(21) 



where the color factor is Cf = (N^ — 1)/2N C — 4/3. Since we neglect the non-leading-log terms, the 
tensor structure of Hfi) q remains the same as H^J in Eq. (0). Similarly, if the final hadron comes 
from the gluon fragmentation, the contribution is 
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rS(l)g 



= W dxtf{x)HiS><x,p, q ) [ 2^ I -C F 1 + {1 Z) D^ h (z h /z) . (22) 

„ J ' JO l T Z7r Jz h z z 

Here z = and £> 9 _>/j(z/j) is the gluon fragmentation function defined as 



dz h 



D g ^ h {z h ) = ^e^ g )d>^ h {z h ,lh) 



= -§=Y, J ^-^ V+/Zh (0\F-»(0Mh)(S,h\F-^y+)\0), (23) 



L h S 

where dg^ h (zh,£h) is given by Eq. (|l2| ) for gluon fields and 

^u(ig)= ^9^>Atg,X) , (24) 

A=l,2 

with s^tg, A) being the polarization vector of a gluon in an axial gauge. 

There are both infrared and collinear divergences in Eq. @ and Eq. (§|). The infrared diver- 
gences come from the phase space where the gluon's fractional momentum goes to zero. These 
divergences will be canceled by the virtual corrections in the quark self-energy diagrams as we will 
discuss later. The collinear divergences when It — > will be absorbed into the renormalized frag- 
mentation functions and are responsible for the evolution of the quark fragmentation function with 
respect to the factorization scale \i. Summing up all the contributions from Eqs. ([l3|), ( ^p| ) and (p2|). 
one has 

dW s f 

' ^ = E / dxtf(x)H$(x,p,q)D q ^ h (z h ,[x 2 ) , (25) 



dz, 

q ~ 

2 f 2 d£ 2 T a s f 1 dz \„ l + z 2 

D q ^ h (z h ,fi ) = D q ^ h {z h ) + 



r 2 di\ 


Q, 


f 1 dz 


10 


2tt 


L h z 



l + (l-z) 2 n 

^F Dg^ h (Z h /Z) 



Cf _ z D q ^ h (z h /z) 
virtual corrections , (26) 



where the renormalized quark fragmentation function D q ^h{zh, M 2 ) satisfies the DGLAP evolution 
equation in Eq. (Q) after inclusion of the virtual corrections. The virtual corrections have to be 
included to ensure unitarity and infrared safety of the final result. We will discuss the virtual 
contributions later in this paper. 



B. Higher Twist Contributions 



In a nuclear medium, the outgoing quark in DIS may experience additional scattering with other 
partons from the nucleus. The additional scattering may induce additional gluon radiation and cause 
the leading quark to lose energy. Such induced gluon radiation will effectively give rise to additional 
terms in the evolution equation leading to modification of the fragmentation functions in a medium. 
These additional terms from multiple scattering are always non-leading-power contributions. There 
are many next-leading- twist contributions to the semi- inclusive hadronic tensor in DIS. One type 
of contribution involves non-leading-twist matrix elements of parton distributions in the nucleus 
which in general are related to two-parton correlations. Since two- parton correlations in a nucleus 
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can involve partons from different nucleons in the nucleus, they are proportional to the size of the 
nucleus and thus are enhanced by a nuclear factor A x l 3 as compared to two-parton correlations in a 
nucleon. In this paper we will consider only those non-leading-twist contributions that are enhanced 
by the nuclear factor. We will neglect those contributions that are not enhanced by the nuclear 
medium. They are in general not related to multiple parton scattering in the nucleus. 




L C R 

FIG. 4. A diagram for quark-gluon rescattering processes with three possible cuts, central(C), left(L) and 
right (R). 



In the LQS generalized factorization scheme for processes in nuclei with large final transverse 
momentum (tfj, ~ Q 2 ), contributions from double parton scattering are leading in otsA 1 / 3 /Q 2 . One 
can similarly apply the generalized factorization to the fragmentation processes in nuclei thanks to 
the LPM interference effect in gluon radiation induced by multiple parton scattering. As we have 
argued in the Introduction and will demonstrate later in this paper, the formation time of gluon 
radiation due to the LPM interference requires the radiated gluon to have a minimum transverse 
momentum i\ ~ Q 2 /MRa ~ Q 2 /A 1 / 3 . The nuclear corrections to the fragmentation function due 
to double parton scattering will then be in the order of a s A x l 3 j ~ a s A 2 / 3 /Q 2 . For large values 
of A and Q 2 , these corrections are leading and yet the requirement l\ <C Q 2 for the logarithmic 
approximation in deriving the modified DGLAP evolution is still valid. Therefore, for large Q 2 
and A, we can neglect processes with more than two parton rescattering and consider only double 
scattering. 

We will work in the LQS framework plfl in which the twist-four parton distributions are factor- 
izable. We will simply apply such factorization of twist-four parton distributions to the study of 
semi-inclusive processes in DIS. In general, the twist-four contributions can be expressed as the con- 
volution of partonic hard parts and four- parton matrix elements jl8| . At the lowest order (processes 
without gluon radiation) in this framework, rescattering with collinear gluons gives the eikonal con- 
tribution to the gauge-invariant leading- twist and lowest-order result in Eq. (|l3|), assuming collinear 
factorization of the quark fragmentation function. As we will show later, rescattering with another 
quark in the leading order will contribute to the semi-inclusive cross section but not to the renor- 
malization equation because there is no collinear divergency. Therefore, we first concentrate on 
processes involving rescattering with gluons. For next-to- leading order processes (with gluon radia- 
tion) involving a secondary scattering with another gluon from the nucleus, shown in Fig. |i| as an 
example, the double scattering contributions to dW^/dzh can be expressed as pll 
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dW, 



D 



dz h 



1 dz 



~ = H / ^-D q ^ h {z h /z) I ^—d Vl dy 



d2 VT ,2, TT D , - 



2tt 



(27T) 



d k T H (y ,y 1 ,y 2 ,k T ,P,q,z); 



x e - i ^*-<A|^(0)7+A+(v 2 -,0 T )A+(tfr,tt r )^(|/-)|i4) . 



(27) 



Here H (y , y 1 , y 2 ikT,P,q,z) is the Fourier transform of the partonic hard part 
H^ u (x, xi, X2, f&r,p, q, z) in momentum space, 



#u„(y ,2/2 ,k T ,p,q,z) = [dx^p-^p-e 



tip II +"2P !/! iBi— a!a)p J/ 2 



2?T 2tt 

*H® v {x,Xi,X2,k T ,p 1 q,z) , 



(28) 



where fcy is the relative transverse momentum carried by the second parton in the double scattering. 
Values of the momentum fractions x,xi, and X2 are fixed by ^-functions and poles in the partonic 
hard part. They normally depend on kr- 

In order to pick up the next-leading-twist contribution, we expand the partonic hard part around 
k T = 0, 



Huv(y ,2/i ,2/2 ,k T ,p,q, z) = H° {y ,y x ,y 2 ,k T = 0,p,q,z) 



1 d 2 H lJV 



fc T =0 



2 dk%dl4 



J\irp ~\~ ... 



(29) 



fc T =0 



This is known as collinear expansion fl22j, On the right-hand-side of Eq. (|29|), the first term gives 
the eikonal contribution to the leading- twist results. It does not correspond to the physical double 
scattering, but simply makes the matrix element in a single scattering gauge invariant. The second 
term for unpolarized initial and final states vanishes after being integrated over kr- The third term 
will give a finite contribution to the double scattering process. Substituting Eq. ( p9| ) into Eq. (p7p, 
and integrating over d 2 kr and d 2 yx 7 we obtain 

" x dz 



dzh 



-D q ^ h (z h /z) 



dy^dy 2 \ (A|^(0) 7 + F+{y 2 ) F+°(y^) ifi q (y~)\A) 



a0 



d 2 =D 



2 dk%dl4 



H uv{y ,2/1,2/2 ,kr,P,q,z) 



(30) 



J k T =0 



where k%A+k%,A + are converted into field strength F a+ F^ + by partial integrations. 



III. QUARK-GLUON DOUBLE SCATTERING 



In this section we calculate the hard part of parton rescattering with gluons. We consider processes 
with a quark and a gluon in the initial state. In these processes, there is first a hard photon-quark 
scattering. The produced parton from the first hard scattering then has a second scattering with 
another initial gluon from the nucleus. We refer to such processes as quark-gluon double scattering 
processes. We defer the calculation of virtual corrections to the latter part of this paper. There are 
a total of 23 cut diagrams corresponding to 3 different double scattering processes, their interfer- 
ences, and the interferences between single scattering and 7 different triple scattering processes (see 
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Appendix). We use the process shown in Fig. |j with three different cuts as an exam/pie to outline 
the steps that are essential to evaluate the hard part of double scattering. 

Since we are interested in the leading log behavior of the fragmentation functions, we will only 
retain the leading contributions in the limit of £t — » (£t being gluons transverse momentum) and 
neglect contributions that are finite when It — 0. This corresponds to the leading log approximation 
for 1% < Q 2 . 



A. Hard vs Soft Rescattering 



in Ref. 18 1, the hard partonic part in Fig. |]with a central cut is 



Following the collinear approximation involving leading twist-4 parton distributions, as outlined 

[| with a 

dxi d,X2 
-i 



H cuv(y ,2/1,2/2 ,kr,p,q,z) 



— ± — ^ c i Xl p + y +ix 2 p + y 1 +i(x-x 1 -x 2 )p + y 2 

2tt 2tt 



d 4 e 



1 

x -Tr 
2 



2wS+{t) 5(1 - z - — ) 



(31) 



Here, 



fr __ Cf_ 4 7 ■ (g + xip) 7 ■ (q + xip - I) 

ap 2N C 9 (q + x lP Y - ze la (q + x lP - If - it lal ' q lp 
x 7- (q + xp-£) 1 -(q + xp) 2irg 

(q + xp — £) 2 + ie (q + xp) 2 + ie q ' 



and 



E a P(l) = -g ali + 



a/3 



n a i + n £ a 



(n ■ If 



(32) 



(33) 



is the polarization tensor of a gluon propagator in an axial gauge (n ■ A = 0) with n = [l,0~,0j_] 
and £, £ q = q + (x\ + x^)p + &t ~ £ are the 4-momenta carried by the gluon and the final quark, 
respectively, z = £~ jq~ is the fraction of longitudinal momentum (the large minus component) 
carried by the final quark. 

To simplify the calculation of the trace part of H ap and extract the leading contribution in the 
limit £t — > and fc^ — > 0, we again use the collinear approximation, 



p a H apP P « 7 • l q 



-Tr 



j-p^H^pf 



(34) 



The two longitudinal components in the ^-integration are fixed by the two ^-functions in Eq. ( J31 
and the remainder becomes d£ T /(Ait) 2 . To carry out the integration over x,X\ and X2, we rewrite 
the ^-function in Eq. (^) as 

1 



s + (e 2 q ) 



2p + q z 



—8(xi + x 2 - x L - x D - xb) 



where xb — Q 2 /2p + q is the Bjorken variable and 



XL 



2p+q z(l — z) 



, xd 



fc| - 2k T ■ £ T 
2p+q~z 



(35) 



(36) 
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With the help of the on-shell condition ( 2 = and £ 2 = in the <5-functions, one can also simplify 
the variables of the propagators in Eq. (fi2|): 

(q + xp) 2 = 2p + q~(x - x B ) ,(q + xp- £) 2 = 2p + q~z(x - x L - x B ) , 
(q + xip) 2 = 2p + q~(xi - x B ) , (q + Xip - I) 2 = 2p + q~z(xi - x L - x B ) ■ (37) 

The integration over x, x\ and x 2 can now be cast in the form 

dx\ dx2 e ixiP+y ~ +iX2P+v i +i ( x ~ Xl ~ X2 ^ p+y 2 



lc{y ,Vi ,V 2 Jr,k T ,p,q,z) = / dx 



2n 2tt [x\ — x B — ie)(xi — xl — x B — ie) 
5(xi + x 2 - x L - x D - x B ) 



(38) 



(x — xl — x B + ie)(x — x B + ie) 

Two of the above integrations can be carried out by contour integration. There are four possible poles 
in the denominators and two of them are used to obtain the residues of the contour integration. There 
are four possible choices for the pair of poles representing subprocesses with different kinematics. 
When we choose the poles 

x = x B + x L , xi = x B + x L , (39) 

which is the momentum fraction carried by the initial quark, the momentum fraction carried by the 
initial gluon in the rescattering is x 2 = Xd, which vanishes when fc^ — > according to the definition 
in Eq. (^). In this case, the final gluon is produced via the final state radiation of the hard photon- 
quark scattering. After the final state radiation, the quark becomes on-shell and encounters another 
scattering with a very soft gluon. This subprocess corresponds to rescattering with a soft gluon after 
a hard collision and is normally referred to as a hard-soft process p8[ . On the other hand, if we 
choose the poles, 

x — x B , xi — x B , (40) 

the momentum fraction carried by the initial gluon in the rescattering is x 2 — xl + xd, which is finite 
(xl) even when kr — > 0. We call this type of rescattering hard and refer the corresponding double 
scattering as double-hard processes. According to such a choice of a pair of poles, the quark becomes 
on-shell immediately after the hard photon-quark scattering. The gluon radiation is actually induced 
by the hard secondary quark-gluon scattering. It is thus produced by the initial state radiation of 
the secondary scattering. The other two combinations of the poles, 

x = x B + x L , Xi = x B ; 
xi = x B + x L , x = x B , (41) 

represent the interferences between hard-soft and double-hard processes. Using the relations imposed 
by the (5-function, x 2 — x B + xl + xd — xi, and the momentum conservation, X3 — x± + x 2 — x, 
one can verify that there is a momentum transfer xl between the initial quark and gluon fields. If 
we consider the quark and gluon come from different nucleons inside the nucleus, this means that 
there is a momentum transfer between different nucleons in these interferences processes. Such an 
observation will become important when we study later the parton matrix elements involved. 
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One can easily find out the residues of the contour integrations with the above four possible choices 
of poles. After making a change of variable x\ + x 2 — xl — xl> — > x, one finds 

f S(x x ) 

Ic{y~,yT,V2^T,kT,p,q,z)= I dx 2 -^-I c (y~,yi,y2jT,k T ,x,p,q,z); 



Ic(y-,Vi, V2,ir, kr, x,p, q, z) = e ^ +x ^ + y~+ iXDP+ ^ -^e(-y^)e(y- - yf) 

x (1 - e - ixtp+y ')(l - e -^iP + fe"-fD) . (42) 

In the above contributions, the interferences between hard-soft and double-hard processes have the 
negative sign relative to the hard-soft and double-hard contributions. Previous studies jl8],^3| have 
considered only hard-soft and double-hard processes, but neglected the interferences. For large It 
or xl , this is justified because the interference terms will vanish due to the oscillation of the phase 
factors in the integration of and y^ over the size of the whole nucleus. Note that in the double- 
hard term the phase factor exp[ixLp + (y^ — y^)} does not vanish after integration over yj~ and y%, 
because color confinement requires that \y± — y%\ remain the size of a nucleon. We will discuss 
other consequences of color confinement later in this paper. At the It — + limit, one can no longer 
neglect the interference terms because they exactly cancel the contributions from the hard-soft and 



double-hard processes. This was pointed out in Ref. |23|. However, our final results, for the first 
time explicitly demonstrate such cancellation. In the latter part of this section, we will discuss the 
physical implications of such cancellation and its connection with the LPM interference effect. It 
will help us to reorganize our final results. 

Because of the collinear (or leading log) approximation, the tensor structure of the double scat- 
tering contributions is generally the same as in the leading order single scattering, 

Kb-.*.*.*.**') = SV.*-.*.*...*.,.). (43) 



-TfD _ _ _ fd£ T a s 1- 

H c {y ,2/i ,2/2 ,kr,x,p,q, z) = J -p- — C F — 



After completing the trace in Eq. ( p4[ ) and using the results of the contour integration in Eq. ([l2|), 
we have the contributions from Fig. |]with a central cut, 

d£' T a s „ 1 + z 2 

27ra s — _ _ _ 
x -jy—Ic{y ,yi,y2^T,kr,x,p,q,z), (44) 

where Ic is given in Eq. © and H$ in Eq. @. 

In addition to the central-cut diagram, we also need to consider asymmetrical-cut diagrams in 
Fig. U that represent interferences between single and triple scattering. The trace part is exactly 
the same as in the central-cut diagram. The differences are that one of the quark propagators is 
replaced by the cut propagator in the central-cut diagram. In addition, each of the asymmetrical- 
cut diagrams has only two possible pairs of poles, whereas there are four choices in the central-cut 
diagram. After similar contour integration, we have for the left-cut diagram 

f dxi dx2 e ixip+v ~ +iX2P+y i +i ( x ~ xl ~ X2 } p+y 2 

lL(y~,yi,y2~jT,k T ,p,q,z)= dx 



2ir 2ir {x\ — xb — ie){xi — xl — xb — ie) 

5(x -x B — x l ) 

(xi + x 2 - x B - x L - x D - ie)(x - xb + ie) 
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= / dx -^ 2 XB ^ I L (y ,y l ,y 2 ,i T ,k T ,x,p,q,z); 



L 

„+„-- 



lL(y-,vT,V;,er,hr,x,p,q,z) = -^^'^^^6^ - % )6{y- - y^) 

-ix L p+ (y~ ~y~) \ 



X (1 



')• 



(45) 



In the second step of the above equation we have made a change of variable x — xl — > x. On the 
left-hand side of the left-cut diagram, the kinematics of the final state gluon radiation in the single 
hard photon-quark scattering requires the quark line to be off-shell, i.e., x = xb- Similarly, on the 
right-hand side, only one of the quark propagators connecting the gluon radiation vertex can be on 
shell. This leads to two possible pair of poles in the propagators, 



xi = x B + x L , x 2 = x D ; 

xi = x B , x 2 — x L + x D , (46) 

giving rise to the two contributions in Eq. (^) . They correspond to soft and hard first quark-gluon 
rescattering, respectively, in the triple scattering process. In both cases, £3 = xd- So the second 
quark-gluon rescattering in the triple scattering process on the right-hand side of the left-cut diagram 
is always soft. 

Similarly, one can also obtain for the right-cut diagram of Fig. ||, 

I R {y-,y^,y 2 ,£ T ,k T ,x,p,q,z) = -e^ +x ^ + y- +ix ^ + ^-y^6(-y^)e(y^ - yf) 



x (1 -e 



-ix L p+y 2 \ 



Other coefficients of the asymmetrical-cut diagrams are exactly the same as in 
gram given in Eq. (0). 



(47) 

the central-cut dia- 



xp 
Ap 











x 3 p + k T 






x 2 p + k T 







2/1 



Hp 



Ap 



FIG. 5. A diagram for the gluon-gluon rescattering process. 



Since QCD is a non-Abelian gauge theory, rescattering can also happen between radiated and 
target gluons. Shown in Fig. ||is an example that involves gluon-gluon rescattering. The contribution 
to the hard part of double scattering from this central cut diagram is 



a s 1 + z 2 



(£ T - k T ) 2 27r 1 - z 
^^7 Fig .|(y" ,y^,y 2 ,£T,k T ,x,p,q,z) 
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I Fls $(y-,yi,y2,eT,k T ,x,p,q,z) = e^ + ^ + y~ +i ^ + ^-y^e{-y^)d{y- - y^) 

which has a structure very similar to the contribution in Eq. (|42] ) and ( ff4|) from the central-cut 
diagram in Fig. |J. The four different terms correspond similarly to the hard-soft, double-hard 
scattering and their interferences. In the hard-soft process, the emitted gluon encounters another 
scattering with a very soft target gluon after it is produced in the final state radiation of the hard 
photon-quark scattering. Note that the initial quark has momentum fraction x — ib+il+id / (1— z) 
while the soft gluon carries momentum fraction X2 = xd — Xd/(1 — z). In the double-hard process, 
the gluon is induced by the secondary hard scattering, similar to the quark-gluon rescattering, but 
in this case via a three-gluon vertex. The kr dependence of the cross section is typical of the 
gluon-gluon interaction. 




FIG. 6. A diagram for the interference between quark-gluon and gluon-gluon rescattering. 

As another example, we illustrate the interference between quark-gluon and gluon-gluon rescat- 
tering processes in Fig. rJ. The contribution from the central-cut diagram is, 



-ttD 



H Fi g Ay >2/i ,2/2 ,k T ,x,p,q,z) = J d£ T 

2na 



2 l T ■ (£t - k T ) a s C A l + z 2 



Nr. 



fc T )2 2tt 2 1-z 
lFi g Wy~ , Vi , V2 , It, k T , x, p, q, z) 



lFi g %(y-,y^,y2,eT,k T ,x, Pl q,z) = -e^+^^^+^tor-^)^-)^- _ y ~) 

x jg^-DP + j/^/(i- 2 ) _ e - ix LP + y2^i — e~ iXLP+ ( y ~ ~ y i^] 



(49) 



Notice that all interference terms in Eqs. ( |45f}47| ) and (j49) have the opposite sign as compared to 
symmetrical diagrams. The calculation of the other diagrams with all possible cuts is tedious but 
similarly straightforward. The splitting functions in all quark-gluon double scattering and their 
interferences have the same form in z as the q — > qg splitting function in the single scattering case. 
The residues of the contour integrations in different processes arc different and the resultant different 
phase factors are essential to give rise to the interesting physical consequences as we discuss below. 
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B. LPM Interference 



Before we list the results of all processes, it is helpful to discuss first the underlying physics in 
gluon radiation induced by quark-gluon double scattering. As we have seen in the calculation of 
the contributions from diagrams in Figs. ||||, the so-called hard-soft processes correspond to the 
case where the gluon radiation is induced by the hard scattering between the virtual photon and an 
initial quark with momentum fraction x b ■ The quark is knocked off-shell by the virtual photon and 
becomes on-shell again after radiating a gluon. Afterwards the quark or the radiated gluon will have 
a secondary scattering with another soft gluon from the nucleus. We denote such a hard-soft process 
by the diagram in Fig. 0, where the off-shell quark is marked by a filled circle. In the double hard 
processes, on the other hand, the quark is on-shell after the first hard scattering with the virtual 
photon. The gluon radiation is then induced by the scattering of the quark with another gluon that 
carries finite momentum fraction xl + xr>- We denote such hard rescatterings by the diagram in 
Fig. D where an off-shell parton is again marked by a filled circle. 




FIG. 7. The diagrammatical representation of hard-soft processes. The off-shell quark line is marked by 
a filled circle. 




FIG. 8. The diagrammatical representation of hard rescattering processes. The off-shell parton lines are 
marked by a filled circle. 



The conventional categorization of soft and hard rescattering by LQS is according to the hardness 
of the virtual partons or the momentum fraction carried by the second initial parton. The hardness of 
the off-shell parton is characterized by the virtuality that is proportional to t\. In the limit £t — > 0, 
these two processes become identical. That is why they cancel each other completely as we see from 
the final results [Eqs. (JI|), (^5|), (f47|), fl48[ ) and j49|)]. Only gluon radiation with finite transverse 
momentum has non-vanishing contribution to the radiation matrix elements. In the collinear limit 
It — > 0, the original terms "soft" and "hard" used by LQS in the large It processes becomes 
ambiguous. But we will keep using them in this paper for historical reasons and for consistency with 
the large £t limit. Physically, one should interpret the soft rescattering as the processes in which the 
gluon radiation is induced by the original first hard scattering (7-quark) , while the hard rescattering 
as the processes in which the gluon radiation is induced by the rescattering, as illustrated in Figs. [7] 



1G 



and The classification of these two different processes according to the pole structure of radiation 
matrix elements is most convenient in discussing the interference effects. 
We define the formation time of the gluon radiation as 

_ 2q-z(l-z) 

T f = p ■ ( 50 ) 

For very collinear (£t — > 0) gluon radiation the formation time can become much larger than the rel- 
ative distance between the two scattering. Then the two radiation processes should have destructive 
interference, leading to the LPM interference effect. One can see that our results [Eqs. ji^), ([l5j), 
(fl7|), ( |48| ) and (Eoj)] have a clear manifest of such an LPM interference effect due to the cancellation 
by the interferences between hard-soft and double hard scattering. Since the gluon's formation time 
should be smaller than the nuclear size MRa/v + ln the given frame, the gluon should then have a 
minimum transverse momentum £ T ~ Q 2 / 'M Ra, as we will demonstrate later in this paper. 

We now reorganize the contributions of different processes according to the above classification of 
hard-soft and double hard rescattering. We will list our complete calculation of quark-gluon double 
scattering in the following. To help understand the reorganization of different processes we provide 
in the Appendix all the radiation amplitudes induced by single, double and triple scattering. There 
we use the approach of helicity amplitudes for high-energy parton scattering and the results are the 
same as our complete calculation in the limit of soft radiation (z g = 1 — z —* 0). 

All the central-cut diagrams are shown in Fig. [| The contributions from double-hard (Fig. ||a) 
and hard-soft scattering (Fig. ||b) are, 

k 2 4. 

x C A =. T =, e" ZXLP {v -»i +y * } , (51) 

e 2 T (e T ~k T ) 2 

e i(x+x L ) P +y-+ix D p+(y-~y-)^^ e ^_ y -^ y - _ y -j 
ix D p+(y--y~+y~)/(l-z) 




C A e T ■ (It - M 



(52) 



a e 2 T (£ T -k T ) 2 j 
respectively. As we have discussed earlier, the gluon radiation in the double-hard scattering is 
induced by the rescattering of an on-shell quark with another hard gluon. Therefore the contribution 
has almost the same form as obtained by Gunion and Bertsch |H](GB) for soft gluon bremsstrahlung 
induced by a single elastic scattering of on-shell quarks, except for the phase factors and the splitting 
function which gives the GB result in the soft radiation limit (z g = 1 — z — * 0). In the processes we 
consider here, the gluon radiation can also be induced by the first hard scattering that produces the 
quark jet. In these so-called hard-soft processes, the final quark or gluon then has a second scattering 
with another soft gluon from the nucleus. In Eq. ([32]), the first term comes from quark rescattering, 
the second term from gluon rescattering and the last two terms come from the interferences between 
the first two processes. 
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o y 2 yi y 

(d) 

FIG. 9. Central-cut diagrams for double-hard 
(a), hard-soft (b) processes and their interfer- 
ences (c) and (d). The diagramatical representa- 
tion of the soft and hard rescattering are shown 
in Fig. § and || 




FIG. 10. Right-cut diagrams that represent 
interferences between single and triple scatter- 
ing. The triple scattering processes involve dou- 
ble hard and soft (a-c) or hard and double soft 
(d-e) processes. The diagramatical representa- 
tion of the soft and hard rescattering are shown 
in Fig. § and || 



There are also two interferences between double hard and hard-soft scattering as shown in Figs, 
and [)|d. Their contributions are, 

x \c A k T - (k T - It) c kr-e T c i XDP +( v -- y7 )/(i-z) \ e -ix LP + y ^ (53) 
\ 2 i\{lT-k T f Gpifr — kx) 2 J 
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H 



D 



OLs l + z2 e i(x+x L )p+y-+ix D p + (y 1 -!),) 2to ^ ( y~)0(y 



x(-). 



2?T 1 - Z 

Ca kr ■ (kr ~ It) 



Ca- 



k T -e 7 



N c 

Jx D p+y 2 /(1-z) 



Vx) 

ix L p + (y~ ~Vi) 



£ T (e T -k T ) 2 "e 2 T (t T -k T ) 2 



(54) 



It is interesting to note that contributions from the double hard processes or gluon radiation in- 
duced by secondary hard scattering and the interferences with hard-soft processes all vanish in the 
collinear limit of the secondary scattering with a gluon (fey — > 0). What remains is the radiation 
spectrum induced by the first hard scattering (photon-quark). As we will show later, combining with 
interferences between single and triple scattering in the same collinear limit, it gives the eikonal con- 



tribution to the next-to-leading order correction of the single scattering [Eq. (21)] corresponding to 
the first term in the collinear expansion in Eq. ( p9[ ) . For finite Ut the above radiation spectra exhibit 
interesting interference patterns which were discussed in detail by Gyulassy, Levai and Vitev j^] and 
Wiedemann || in the framework of the GW static color-screened potential model. In the present 
paper, we want to study their contribution to the higher- twist corrections to the fragmentation 
functions. 

To complete our calculation we also have to consider all the interferences between single and 
triple scattering. These correspond to the asymmetrical-cut diagrams. The right-cut and left-cut 
dialgrams in Fig. ^ are just two examples. After similar reorganization of all different processes, 
we list in Figs, [h] all the possible right-cut interference diagrams. The contributions from the first 
three right-cut diagrams in Fig. |o[ which involve double hard scattering, are 



RH 



= / Mr 



a, 1 



z 



i(x+x L )p y +ix D p {y 1 -y 2 ) 



x (-) 



2vr 1 



z 

■Ct) 



N, 



k T f 



C A 



,i(x a D -x D -x L )p+ [y 1 -y 2 ) 



0(-%)%2~ -vl) 

Ca 
2 



C A e 



-i(l-z/(l-z))x D p + (y 1 -y 2 ) 



-ixLP + y 2 



(55) 



where x° D = k^/2p + q^ . We have made the variable change kr — > —kT in some of the contributions 
in order to obtain a more compact form of the final result. The contributions from the two hard-soft 
processes in the right-cut diagrams in Fig. mi are 



H RS = 



a s 1 



a i(x+x L )p + y +ix D p + (y 1 -y 2 ) 



2tt 1 

l r 



N c 



C F + C A e 



-i(l-z/(l-z))x D p+(y 1 -y 2 ) 



l T ■ {l T - k T ) C A 



,ix D p+y 2 /(l-z) 



I _|_ e — i ( 1 ~ z /( 1 - z ))xDP + (y 1 -y 2 ) 



(56) 



k T ) 2 2 , 
Similarly, contributions from the left-cut diagrams, which are just the complex conjugates of the 
right-cut diagrams in Fig.nO, are, 



x (-) 



2tt 1 
1 ■ (k T 



z 
■It) 



N r . 



£ 2 T (£ T - kr) 
(j Ae -i( l ~ z /( 1 - z ))xDP + (y^-y 2 ) 



Ca l (x° D -xo-x L )p+(y-- v -) _ Ca 
2 "2 



Q(y - Vi W(vi - y 2 ) 



-ix L p+(y -j/j ) 



(57) 
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H D = f d£ 2 as 1 + Z c i(x+XL)p+v ' +iXDP+(v7 ^ ) 2 ' Kas 
LS / t 2tt l-z N r 



0(y - v\ W(vi - V2 ) 



x (-) 



£ T ■ (£ T - k T ) Ca cIXdp+(v -^ v7)/(1 _ z) 
i 2 T (f T -k T ) 2 2 



I _|_ e -i(l-*/( 1 -*)) a: £>p + (y 1 -y 2 ) 



(58) 



C. Collinear Expansion 

Following the procedure of extracting the next-leading-twist contributions to the semi-inclusive 



DIS cross section, as outlined in Sec. [IB, we now expand the hard part in kr according to Eq. (|29|). 
To simplify our notation, we factor out all the 9- functions and /cy-independent part of H c R L (the 
sum of all the contributions from central-cut, right-cut or left-cut diagrams) and define 

~H D = [ d£ 2 — 1 + z2 c i( - x + x ^r> + y-+ixDP + (v7-y^) 27ras 
' J T 2ir l-z N c 

x [Hge(-y^)9(y- - yf) - H°d(-y^)6(y; - yf) - fffflfcr - y^)9(y^ - %)] . (59) 



From Eqs. (|51|)-(|5q), we have 

H D c {k T = 0) = H D R {k T = 0) = H°( kT = 0) = ^f. (60) 

They all come from hard-soft processes. As we have pointed out before, double-hard scattering 
processes and their interferences all vanish when kr = 0. Note that the combination of ^-functions, 

dyidyz - Vi) + 0{y- - y^)6{y^ - - 6{-y^)6{y- - )] 

dy ± / dy 2 , (61) 



'o Jo 

is an ordered integral limited by the value of y~ . We will refer to any terms that are proportional to 
the above combination of ^-functions as contact contributions. Combining with Eqs. (|27]) and (p3|), 
the first term in the kx expansion in Eq. (|5^), H D {kr = 0), gives the eikonal contribution to the 
next- leading-order correction of single scattering in Eq. ( pl| ) . Such eikonal contributions involving 
non-physical gauge fields do not correspond to any physical double scattering. They only make the 
final results gauge invariant and can be gauged away. However, we do see the importance of including 
all the interferences (right-cut and left-cut diagrams) between single and triple scattering. This is 
essentially a special case of the generalized proof of factorization of leading-twist contributions in 
DIS pb]| where contributions from any number of soft gluon rescatterings can be eikonalized. 

The dominant contributions to the double quark-gluon scattering come from the quadratic term 
in the kr expansion of H D in Eq. (^). We note that these high- twist terms generally involve 
two additional spatial integrations with respect to y-j" and y 2 ■ If no restriction is imposed, these 
integrations will only be limited by the nuclear size and thus will produce nuclear enhancement 
as compared to DIS with a nucleon target. Because of the rapid oscillation of e ixp y , any term 
proportional to Eq. ( |6l| ) that limits the integration over y^ and y 2 to the value of y~ will not 
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have nuclear enhancement. Neglecting any term that is proportional to Eq. (|6l|) (contact term) and 
keeping the leading term when £t — ► 0, we have 



V 2 kT Hg\ kT = = ^(1 - e --^)(l - e -W(if--»D) + 0(x B /Q*t%.) , 
VL fff | fcT =o = + 0{x B /QH 2 T ) , 



V 2 fer ^| fcr=0 = + 0{x B /QH 2 T ) . (62) 

The four terms in correspond to hard-soft, double hard scattering and their interferences. 
Substituting the above in Eq. (E§) and (^), we have the leading contribution to the semi- inclusive 
tensor from double quark-gluon scattering with quark fragmentation 

E^=J2j dxHjg>(xp,q) J^D q ^ h (z h /z)^-C A ^ 
d£l 2ira 



J ^1 2 ^ S j gg( a: ' XL ) + ( virtual correction) , (63) 



where 



i (A|^(0) 7 + F+fo) F+°{y^) ^ q (y-)\A)0(-y^)9(y- - yf) (64) 



is the quark-gluon correlation function which essentially contains four independent twist-four parton 
matrix elements in a nucleus. Since xl as defined in Eq. ( |36| ) depends on £t and z, the matrix element 
T^(x, xl) has an implicit dependence on z and £t< The contribution from gluon fragmentation is 
similarly 



dz h 



2 



x f ^^ ^ 7T ™ s T i £ ) (x,xl) + (virtual correction) . (65) 

The above contribution from double quark-gluon scattering is very similar to the next-to-leading 
order single scattering in Eq. (|2l|). Even the splitting function has the same form except for the 
color factor. The contribution is, however, proportional to the twist-four parton matrix elements. 
Since there is not much restriction on the spatial integral, such twist-four matrix elements will give 
the nuclear enhancement we are looking for. We will return later for more discussions on the A 
dependence of the double scattering contribution. 

In the process of a collinear expansion, we have kept £t finite and took the limit kx — > 0. As 
a consequence, the gluon field in one of the twist-four parton matrix elements in Eq. (|64| ) carries 
zero momentum, corresponding to the soft-hard process. Later in this paper we will discuss how to 
estimate these twist-four matrix elements. Assuming a factorized form, they will be proportional 
to the product of quark and gluon distributions in a nucleon. However, the gluon distribution 
xf g (x) at x — is not defined in QCD. This is because we dropped higher order terms of the 
collinear expansion in Eq. (|6^). These terms are considered higher-twist contributions in the LQS 
framework. As a remedy to the problem, we simply resum a subset of the higher-twist terms in 
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the collinear expansion and restore the phase factors such as exjp(ixTP + y~) which is related to 
the intrinsic transverse momentum of the initial partons, xt = (k^) / 2p + z . In this case, soft 
gluon fields in the parton matrix elements will carry a fractional momentum xt ■ Similarly, one can 
also restore the collinear structure of the scattering amplitude by replacing 1/4 in Eq. ( j62| ) with 

i/4(4 + <*£»• 




L R 

FIG. 11. Diagrams for the virtual correction to single scattering with two possible cuts. 



D. Virtual Corrections 



So far we have not considered virtual corrections which will ensure the final result to be infrared 
safe. The calculation of the virtual correction to the single scattering in Fig. [fl] is very similar to 
the real correction in Fig. |^. One can easily find, 



dz h 



= -£/ dxtf(x)H<?J{x,p >q )D^ h (z h ) £ dzC F \ 



1 + z 2 



(66) 



One can see that the integral over z is infrared divergent at z = 1 or z g = 0. However, such 
divergency cancels exactly the infrared divergency in Eq. (|2l]) of the radiative contribution to the 
single scattering process. The sum is then infrared safe. Using the definition of the +functions [M, 



(l-z). 

with F{z) being any function which is sufficiently smooth at z = 1, one can rewrite 
radiative and virtual correction [Eq. ww] in the following form 



(67) 
the sum of 



df 



D q ^ h (z h , fj?) = D q ^ h (z h ) + 

+ j q ^ gq (z)D g ^ h (z h /z)\ 

The splitting functions are defined as 



1 dz 



-IT^ZI — hq^qg(z)Dq^h{Zh/ 



Iq^gq(z) = Iq^qg (1 - z) 



— + -(5 1 - z) 

(l-z)+ 2 



(68) 

(69) 
(70) 
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Therefore, with the definition of the +function, the (^-function terms in the splitting functions 
take into account the self-energy virtual correction which cancels the infrared divergences from the 
radiative processes. The final renormalized quark fragmentation D q ^h{zh, /J 2 ) satisfies the DGLAP 
evolution equations in Eq. ([[]). 

When cast into the DGLAP evolution equation in Eq. (|l]), the splitting function from the real 
correction can be interpreted as the probability for the quark to radiate a gluon with momentum 
fraction 1 — z. Then one must also take into account the probability of no gluon radiation in the 
evolution to ensure unitarity. Such unitarity requirement gives rise to the same virtual correction as 
calculated from the diagram in Fig. [ll]. In the double scattering case, we will use the same unitarity 
requirement to obtain virtual corrections. The virtual contribution to the quark fragmentation in 
double scattering processes is, for example, 



W$ vU 
dz h 



= -J2 I dxH$(xp,q)D q ^ h (z h )^-C 



dz 



1 + z 2 



<2 ™ S T>,^). (71) 



Nr. 



One can single out the infrared divergent part by rewriting the integral 



1 + z 

dz ~, T aq(. x i x L) =T 'Jx,X L )\ x= l 



dz- 



l-z 



XT A Jx,l 2 T ) 



AT A (x,e 2 T ) = [ dz-?— [2T A (x,x L )\ z=1 - (1 + z 2 )T A {x,x L )) 
Jo 1 — z 



(72) 



The second term is finite since T A (x, xl) is a smooth function of z. The first term can be combined 



with the radiative contribution in Eq. (|63j) to cancel the infrared divergency. With the help of the 
+function, the final result can be expressed as 



W Dq 
dzh 



E J dxH^P^ J j£ J zh U iD q ^( Zh /z) 



d( 2 , 



X 2^ CA 



{\-z\ 



T g A Jx,x L )+S(z-l)AT g A g (x,£ 2 T ) 



(73) 



Here the implicit z-dependence of T a (x,xl) plays an important role in the final result. The above 



integrand will be proportional to the splitting function for single scattering in Eq. ( p9\ ) if one ignores 
the z dependence of T A (x, Xl). Similarly, the final result for contributions from gluon fragmentation 



W D -9 

dzh 



E I dxHjPJ { xp, q) 2 -^ r-^f I ^ h (z h /z) 



N c 



dl 2 



dz 



1 + — T A (x, x L ) + S(z)AT A (x, 4) 



x ^Z Ca 

Z7T 



where we have used the fact that xl in Eq. ( |36|) is invariant under the transform z 

r A ( 



is T A (x,x L ). 



(74) 
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FIG. 12. The diagrams for leading order quark-quark double scattering 
IV. QUARK-QUARK DOUBLE SCATTERING 



We have so far only considered quark-gluon double scattering in a nucleus. After the first photon- 
quark hard scattering, the leading quark can also rescatter with another quark from the nucleus. 
Such quark-quark double scattering processes also contribute to the semi- inclusive DIS at twist- four. 
Unlike quark-gluon rescattering, quark-quark double scattering even contributes at the lowest order 
without gluon radiation. At the lowest order there is only one kind of quark-quark double scattering 
diagram, as shown in Fig. [l2] and its crossing variations. One can easily calculate their contributions 
and obtain 



dzh 



1 



■*C F ^H$(x,p,q) {D g ^ h (z h )6(-y^)d(y- - 



D q ^ h {z h ) [9(y- - y^)6(y- - y-) + 9(-y^)0(y^ - yf)] } , 



(75) 



where to combine all the crossing diagram we use the fact that fields on the light-cone commute 
with each other J2]J| . In the above contributions, the term proportional to the gluon fragmentation 
function is from the central-cut diagrams while the terms containing quark fragmentation function 
are from left and right-cut diagrams (or interferences). Using Eq. ( |6l|) and neglecting the contact 
term, we have 



dW, 



D(0) 



dz h 



E 

9 



dxT^Xx 



2na 6 



8Cj 



0* 



H$(x,p,q) [D g ^ h {z h ) - D q ^ h {z h )] 



where 



^^r^e" p+ ^(A|^(o)^v- ? (^)^(yr)^^fe)|A>e(-^)0(y- 



(76) 



Vi) (77) 



is a four-quark matrix element in a nucleus. This twist-four contribution is explicitly suppressed by 
l/Q 2 as compared to the single scattering case. Such leading-order and high- twist contributions from 
quark-quark double scattering essentially mix the quark and gluon fragmentation functions. In this 
case, there is no induced radiation and thus no energy loss. However, these contributions will change 
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the final differential cross section or semi-inclusive spectrum, since quark and gluon fragmentation 
functions are different. One thus can consider the case as modification of fragmentation functions 
without energy loss. Furthermore, these processes do not change the integrated total cross section. 
One can verify this by using the momentum sum rule J2h I dz hZhD q ^ g ^h{zh) = 1- This is consistent 
with a general theorem that final state interaction will not change the total cross section. 




FIG. 13. An example diagram for next-to-leading order quark-quark double scattering. 



We also note that contributions from quark-quark double scattering are free of any divergences, 
especially collinear divergency. Therefore, they will not contribute to the QCD evolution of the 
effective parton fragmentation functions in the nuclear medium. 

To the next order we have to consider radiative corrections to the processes in Fig. [l^. In addition 
we also have to consider the diagram in Fig. [Tj| Again, calculations of these diagrams are tedious 
but straightforward. Take the diagram in Fig. [l3]for example: The pole structure is exactly the same 
as the central-cut diagram of quark-gluon double scattering in Fig. |J. The resultant contribution is 

dW Fig ] 13^ _ST f 1 dz r> ( I \ [ M T a s n 1 + X B 



dZh ~^J Zh z D ^ {zh/z) J ^ CF \i- z y 



9 



dxT^(x,x L )H^(x,p,q) , (78) 



where 



T q f I \x,x L ) = Y / J ^-dy^dy^e^ +x ^ + y'(l-e-^^ )(l 



e ~ix L p + (y 



x (AHM^My^ivT^^MMei-y^eiy- - y r) ■ (79) 

The four terms with different phase factors in the above equation correspond to hard-soft, double 
hard and the interferences. Notice that the two diagonal terms arc defined exactly the same as 
T^q I \x). This structure is exactly the same as in the quark-gluon double scattering. However, this 
contribution is only proportional to versus in the quark-gluon double scattering. It is thus 
suppressed by l/Q 2 - This remains to be the case for all quark-quark double scattering processes. 
Since we are interested only in the collinear behavior of the gluon radiation processes in order to 
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study the QCD evolution, we will neglect in this paper all radiative contributions from quark-quark 
double scattering. Similarly, contributions proportional to l/^fr from gluon-gluon double scattering 
can also be neglected to this approximation, e.g., in the collinear expansion in Eq. (62). 



V. MODIFIED QUARK FRAGMENTATION FUNCTION 



For a complete result of the semi-inclusive cross section of DIS off a nucleus, one should also 
include the higher-twist contribution to the quark distributions and their QCD evolution as studied 
by Mueller and Qiu (^6). Here we will simply replace the leading-twist quark distributions f£(x) by 
fq{x,ix 2 j) which contains nuclear modification to the quark distributions and their QCD evolution 
and /ij is the factorization scale for the quark distributions in a nucleus. We should note that even 
the leading- twist quark distributions f q (x) in a nucleus are different from A free nucleons [f q {x)\. 



A. Modified Evolution Equations 



Including the higher-twist contributions to the quark distributions and summing up all the leading 



contributions from single and double scattering processes in Eqs. fl6q), ( |73| ) and (74), we have 
dW MV = 



dzh 



]T J dxtf(x,tf)H$(x,p,q)D q ^ h (z h ,ix 2 ) 



(80) 



as the total semi-inclusive tensor in DIS off a nucleus up to twist- four corrections. We define the 
modified effective quark fragmentation function as 



~ d£ 2 a f 1 dz 

D q ^ h (zh, M 2 ) = Dq^h(zh, A* 2 ) +/ -J-TT —[^ q ^ qg (z,x,x L ,£^)D q ^ h (z h /z) 

Jo l T ^ Jzh Z 

+ Aj q ^ gq (z,x,x L ,eT)Dg^ h (z h /z)] , 
where D q _,f l (zf l , /j 2 ) is given in Eq. (pojj for leading-twist contributions and 



A~f q ^ qg (z,x,x L ,£T) 



1 + z 2 



. . )+ Tf g (x, x L ) + 6(1 - z)AT q A g (x, e T ) 



{t% + {k 2 T ))N c mx,fi 



(81) 



(82) 
(83) 



|) and AT^(i,4) is given in Eq. ©. 



Aj q ^ gq (z,x,x L ,£^) = Aj q ^ qg (l - z,x,x L ,£ 2 T ). 

The twist- four matrix element T q (x, xl) is given in Eq. 
Since the two-parton correlation function T^ g can involve partons from two different nucleons inside 
the nucleus, it should have a nuclear enhancement that is proportional to the nuclear size Ra, or 
Tqgl fq x A 1 /^ , as we will show in the next subsection. Because of the LPM interference structure 
contained in T q A g (x,XL) — (xl ~ £\ ~~ * 0); om y those radiated gluons induced by multiple parton 
scattering that have finite transverse momentum, l\ ~ Q 2 /A x ^ 3 contribute to the fragmentation 



function. Therefore, the nuclear correction to the fragmentation function in Eq. 81 is proportional 
to a s A 2 / 3 /Q 2 from double parton scattering. For large A and Q 2 , this is the leading nuclear 
correction. Since T qg (x, Xl) / f q {x) is proprotional to the gluon density per unit transverse area in 
the nucleus p g and a s /Q 2 can be considered as the parton rescattering cross section a g , the leading 
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nuclear correction to the fragmentation function here is then proportional to a g p g R,A- This is exactly 
the same as the expansion parameter in the opacity expansion approach in Refs. 

We should emphasize here that the factorized form of the semi-inclusive tensor in Eq.(|80|) only 
serves to define the effective quark fragmentation function D q ^h(zh, fi 2 )- Such a modified fragmen- 
tation function for final state particle production has an explicit dependence on the initial parton 
distribution through the high-twist double scattering processes. Therefore, the factorization for 
semi-inclusive processes in DIS is broken explicitly at twist-four correction. This is a natural conse- 
quence of the non- vanishing parton energy loss at twist-four when the leading quark suffers multiple 
scattering through the nuclear medium. 

Taking the derivative with respect to the collinear factorization scale /i 2 , we obtain the modified 
DGLAP evolution equation in leading order of a s for the modified quark fragmentation function, 

rl dz r 



dD q ^ h (z h ,fi 2 
din ji 2 



Oh 

2tt 



%-> qg {z,x,x L ,ii )D q ^ h (z h /z, n ) 



+ %^ gq (z, x, x L ,n 2 )Dg^ h {zh/z, is 2 )] 



(84) 



The modified splitting functions are defined as 



~1 q ^ qg {z,x,x L ,ii 2 ) = j q ^ qg (z) + A~f q ^ qg (z,x,x L ,ii 2 ) 

Jg^ gq (z,X,X L ,H 2 ) = 7,->g fl (l ~ Z, X , X L , ^ ) , 



(85) 
(86) 



where ^f q ^ qg {z) given in Eq. ( p9[ ) is the splitting function in single scattering processes. If one 
only considers single-jet events in DIS, the gluon fragmentation function does not contribute in the 
leading order. Before we carry out similar calculations for the higher-twist correction to the DGLAP 
evolution of the gluon fragmentation functions, we may assume that the gluon fragmentation function 
which enters into the above equation follows the normal DGLAP evolution, 



dD g ^ h (z h ,iJ 2 ) 
d In /i 2 



OLs_ 

2tt 



1 dz 

, z 



7g^qq( z ) D q^h{zh/ z, V' 

9=1 



) + 7 3 -99 (z)D g ^ h {z h /z,H 2 ) 



(87) 



where the normal splitting functions are 

79-99-= ^[* 2 + (l- Z f]> 
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J -n f )S(z-l) 



(89) 



and nf is the number of quark flavors. All the fragmentation functions obey the momentum sum 
rule, 



dz^z Da^^z^n 2 ) = 1 

h 



(90) 



One can check that the modified quark fragmentation function D q ^h(z, p 2 ) in Eq. ( [si] ) still satisfies 
the momentum sum rule. Such a momentum sum rule might seem count-intuitive since there is 
momentum transfer of x lP + to the fragmentating quark from the second partons in the nucleus in 
the double-hard processes. However, one should note that the initial quark from the nucleus carries 
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momentum xsp + in this case. In single scattering and hard-soft processes, however, the initial 
quark carries momentum (xb + xl)p + - So the total momentum transfer from the nucleus to the 
quark is the same in both single and double scattering. The momentum sum rule for the modified 
fragmentation functions should still be valid. 

Compared to Eq. (Q) , the renormalization equation for the modified quark fragmentation function 
in Eq. ( p^ ) is similar to the original DGLAP evolution equation for the fragmentation functions in 
vacuum. However, the modified splitting functions 7 [Eq. ( p5| )] have an extra term A7 from induced 
gluon radiation. Similarly to the ordinary parton cascade in vacuum, the induced radiation will 
soften the modified quark fragmentation function. As we have argued in the Introduction, such 
softening will be the only experimentally measurable consequence of quark energy loss in a medium. 
Another important difference between the modified and the original DGLAP evolution equations is 
that the induced splitting function A7 depends on the twist-four parton matrix elements T^(x, xl) 
of the nucleus. Because of the interferences between hard-soft and double-hard processes, T^(x, xl) 
explicitly manifests the LPM effect which modifies the gluon radiation spectra such that the induced 
splitting functions are also modified from their form in vacuum. The induced splitting functions also 
explicitly depend on the factorization scale fi. 

B. Twist- four Parton Matrix Elements 

As we have observed, modified DGLAP evolution equations for parton fragmentation functions 
depend on both the parton distribution f A (x,nf) and two-parton correlation functions T^(x,xl)- 
Naturally, they are essential to any numerical solution to the modified DGLAP evolution equations 
for the fragmentation functions. One also has to know them, especially their dependence on the 
nuclear size, in order to have any numerical estimate of the effective quark energy loss. 

As defined in Eq. (|64|), the parton correlation function T^(x : xl) essentially contains four inde- 
pendent twist-four parton matrix elements in a nucleus. They arc in principle not calculable and 
can only be measured independently in experiments just like parton distributions. However, under 
certain assumptions, one might be able to relate the parton correlation functions to single parton 
distributions in the nuclei and in the meantime obtain the A-dependence of the parton correlations. 
There are four parton field operators and three spatial (in the longitudinal direction) integrations 
in T^(x, xl)- We assume that the nuclear wave function can be expressed as a multiple- nucleon 
state and each nucleon is in a color singlet state. In this case, there should not be long range color 
correlation between two nucleons due to color confinement. In other words, the interaction between 
nucleons inside a nucleus can only be mediated via color-singlet objects. Consider first the case 
where the quark and gluon fields operate on different nucleons inside the nucleus. This means that 
the double scattering happens with two partons from two different nucleons. Because of the color 
confinement, y\ and j/2 in Eq. ( |64| ) are restricted such that \y\ — 2/2 j < r/v, where rpj is the radius 
of a nucleon. The integrations over y\ and j/2 should, therefore, give the length scale of t^Ra, 
where Ra — 1.12 A 1 / 3 is the radius of the nucleus. In this case the correlation functions in Eq. ( |64| ) 
should be approximately proportional to A 4 / 3 . If both the quark and gluon fields operate on the 
same nucleon state inside the nucleus, y\ and j/2 are both restricted to be within the size of the 
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nucleoli. Then the correlation functions should only be roughly proportional to A. This corresponds 
to double scattering with two partons in the same nucleon. As we have discussed before, we will 
neglect these contributions because they yield no nuclear enhancement. 

Among the four independent matrix elements in T^(x,xl), there are two diagonal terms, one 
proportional to a phase factor exp[i(a; + XL)p + y~ + ixTP + (Vi ~ 2/2")] an< ^ another to exp[ixp + y~ + 
i{x,h + xt)p + {Ui ~ 1/2)]- The first term corresponds to hard-soft double scattering processes in 
which the gluon carrying momentum fraction xt is soft. The second term comes from double- hard 
scattering in which the gluon carries large momentum fraction xl +%t- In this paper, we adopt the 
approximation used in Ref. Jlq ). Under such an approximation, one simply relates the two-parton 
correlations to the product of two single parton distributions, 

j..— j..— „ixip + y~ +ix 2 p + {y± — j/J ) 



. dy l dy 2 



x \{A\^ q {Q) 1 + F+{yz)F + %y^)4 !q (y-)\A)6{-y 2 )0{y--y^) 



C_ 

xa 



^f A (x 1 )x 2 f?(x 2 ). / (91) 



where xa — 1/MRa, fq{%) is the quark distribution inside a nucleus as defined in Eq. (14), and 
fg (x) is the gluon distribution inside a nucleon as defined by 

${x) = / ^-e^ + y-(N\F a + (0)F+°(y-)\N) . (92) 

C is assumed to be a constant, reflecting the strength of two-parton correlation inside a nucleus. As 
we have argued, C should be proportional to the nucleon radius rjv, characterizing the confinement 
scale. We also assume that the nuclear distribution takes a Gaussian form p{r) ~ exp(— r 2 /2R 2 A ). 

2 _2 _ 

In terms of the light-cone coordinates, we have then p(y~) = poexp(y~ /2R A ) where R A = 
V2RaM/p + and M is the nucleon mass. The distribution is normalized such that 

_ 2 -2 

dy-p e- v I 2R a = 1. (93) 

With such a distribution, the third integration in Eq. ( |9l| ) gives J dy 2 ~ R A ~ l/x A p + - 

In addition to the two diagonal terms, there are also two off-diagonal matrix elements in T^ g (x, xl) 
which come from the interferences between hard-soft and double-hard scattering processes. These 
two new matrix elements introduced by the interferences have never been studied before in DIS. 
Similar matrix elements have been discussed before in the context of pA collisions ]27| ]. For finite 
value of xl, this involves off-diagonal parton matrix elements and they should be related to the 
so-called skewed parton distributions of nucleons |28| . For a rough estimate here, we can generalize 
the approximation on parton correlations in Ref. Jl8| to the off-diagonal matrix elements. As we 
have pointed out in the calculation of the cut-diagrams in Fig. ^, in the interferences between double- 
hard and hard-soft processes, there is actually momentum flow of xlP + between the two nucleons 
where the initial quark and gluon come from. Without strong long range two-nucleon correlation 
inside a nucleus, the amount of momentum flow xlP + should then be restricted by uncertainty 
principle. Note that the phase factors in the off-diagonal matrix elements have an additional factor 
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of exjp(±ixLP + V2 ) relative to their diagonal counterpart. Assuming the same Gaussian nuclear 
distribution in light-cone coordinates and using 

dy- p e- y ~ 2 l 2R ~* MxhP+y ~ = e' x ^ /x ^ , (94) 



we should have a similar approximation for the off-diagonal matrix elements: 



T q A g (x,x L ) = —lff(x + x L ) x T f g N (x T ) + tf(x)(x L + x T )f»(x L + x T )](l - e-*d*A) . (96) 



^ dn~d.tr P ix ^P + V +i*2P + {y 1 -y 2 )±ix L p + y 2 

x l -{A\%{Q) 1 +F+{y-)F+%y-)%{y-)\A)e{-y-)e{y--y-) 

= —tf{ Xl )x 2 f»{x 2 )e-^/ x -. (95) 

Relative to the diagonal matrix elements, the off-diagonal ones are suppressed by a factor 
exp(— x 2 L /x 2 A ). Combining all the four terms together, we have 
C_ 

XA ' 

Notice that ry = \/xlp + is the gluon's formation time. Thus, xl/xa = ^ j ~aI t } with — RaM/p + 
being the nuclear size in the our chosen frame. It is then clear from the above that the interferences 
between double-hard and hard-soft scattering cancel exactly the double scattering contributions for 
collincar gluon radiation whose formation time is much larger than the nuclear size. The LPM effect 
now explicitly manifests itself via the effective matrix elements in the gluon radiation spectrum that 
is induced by double parton scattering. 

The diagonal parton matrix elements that arise from hard-soft and double hard processes are 
exactly the same as what appear in processes with large final transverse momentum such as di- 
jet production studied in Ref. |Ts|| . These diagonal matrix elements have an enhancement linear 
in nuclear size Ra relative to the single parton distribution in nuclei as in Eq. (|9l]). The off- 
diagonal parton matrix elements, on the other hand, come from the interference between hard-soft 
and double hard processes. Because such parton matrix elements involve transferring momentum 
fraction xl between two nucleons inside a nucleus, they should be suppressed when xl is larger 
than what the uncertainty principle allows. For a Gaussian nuclear distribution, these off-diagonal 
matrix elements are exponentially suppressed for large xl/xa as in Eq. (|95|). This is why one can 
neglect the interferences between soft and hard rescattering in the LQS Q study of processes with 
large final transverse momentum, It » Q 2 /A 1 ^ 3 . However, in our current study of fragmentation 
processes, the off-diagonal matrix elements in the interferences are very important. The interference 
terms cancel exactly the normal contributions that contain the diagonal matrix elements as shown 
in Eq. ( |96|) for £t = 0. For small but finite £t, the cancellation is not complete. The effective parton 
matrix element (diagonal minus the off-diagonal ones) is then proportional to x 2 L /x\ which depends 
cubically on the nuclear size Ra- It, however, also vanishes quadratically with xl ~ &[<■ Therefore, 
the interferences effectively cut off the transverse momentum of the radiated gluon at i\ ~ Q 2 /A 1 / 3 . 
In terms of the LPM interference, this means that the formation time of induced gluon radiation 
cannot be larger than the nuclear size Ra- Such a consequence of the LPM interference effect is 
critical for us to apply the LQS generalized factorization to the problem of nuclear modification of 
fragmentation function in this study. 
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C. Parton Energy Loss 



So far we have derived the modification to the quark fragmentation functions and their DGLAP 
evolution equations. Given the twist-four parton matrix elements, one can then solve the modified 
DGLAP evolution equations for the modified fragmentation functions. Doing so, one can effectively 
resum all the leading-log corrections at twist-four. We leave the numerical solutions to future studies. 
In this paper, we instead will formulate the effective energy loss of leading quarks. 

In principle the modification of the fragmentation functions would be the only experimental effect 
of induced gluon radiation via multiple scattering. One can never directly measure the energy loss 
of the leading quark. The net effect of the energy loss is the suppression of leading particles on 
one hand and the enhancement of soft particles on the other, leading to the modification of the 
fragmentation functions. One can then experimentally characterize the parton energy loss via the 
momentum transfer from large to small momentum regions of the fragmentation functions. We will 
discuss this in detail when we numerically calculate the modified fragmentation functions. 



Upon a close examination of Eq. (81), we see that the first term is the renormalized fragmenta- 
tion function in vacuum. The rest are particle production induced by the rescattering of the quark 
through the nuclear medium. In particular, the last term is particle production from the fragmen- 
tation of the gluon which is induced by the secondary scattering. Such particle production is at the 
expense of the energy loss of the leading quark. We can thus quantify the quark energy loss by the 
momentum fraction carried by the induced gluon, 

^ mi r 1 



(Az g )(x B ,^ 2 ) = J -jr J dz^-z Aj q ^ gq (z,XB,x L ,£?r) 



" t dz\l + (1 - z?] T ^b,x l ) C A al 

4 Jo 1+1 n 7#**,M?)*c(4 + <**»' 

Using the approximation of the parton correlation functions in Eq. ( |96| ) and changing the integration 
variable, we have 

9 N c f q A (x B ,^ 2 )x A Q 2 Jo z(l-z) J x L (x L +x T ) y 

x [fq (x B + x L ) xtJ^xt) + f£(x B ){x L + x T )f"{x L + x T )) , (98) 

where x^ = n 2 /2p + q~ z(l — z) = x B /z(l — z) if we choose the factorization scale as fi 2 = Q 2 . 

One can numerically evaluate the above in the future with assumed parton distributions in nuclei. 
It is nevertheless interesting to discuss the nuclear dependence of the parton energy. The exponential 
factor in the above equation comes from the LPM interferences. It regularizes the integration over 
xl which is otherwise divergent. Combining with the QCD radiation spectrum, it also limits the 
integration over xl to xl < xa and gives J dxLjx\ ~ 1/xa- The fractional energy loss by the quark 
is, therefore, proportional to 

Since xa = 1/MRa., the energy loss thus depends quadratically on the nuclear size. The extra size 
dependence comes from the combination of the QCD radiation spectrum and the modification of 
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the available phase space in It or i/, due to the LPM interferences. One also notes that, though 
the fractional energy loss is suppressed by 1/Q 2 for a fixed value of xb, the total energy loss, 
AE = q~(Az g ), is not. 

VI. SUMMARY AND DISCUSSIONS 

Working in the framework of generalized factorization of higher- twist parton distributions, we have 
studied multiple parton scattering in deeply inelastic eA collisions, in particular the induced gluon 
radiation and the resultant modification to the final hadron spectra. We have defined modified 
quark fragmentation functions to take into account the effect of multiple parton scattering. We 
have presented a detailed derivation of the modified quark fragmentation functions and their QCD 
evolution equations to the next-leading- twist. The modification is shown to depend on twist-four 
parton matrix elements, both diagonal and off-diagonal, in nuclei. 

Depending on whether the gluon radiation is induced by the secondary scattering, one can cate- 
gorize the multiple parton scattering as soft or hard, according to the fractional momentum carried 
by the secondary parton involved. We have considered both soft and hard scattering and their 
interferences. We have shown that these are exactly the so-called LPM interference effect. We have 
demonstrated that LPM interferences modify the available phase space in the emitted gluon's mo- 
mentum. Coupled with the gluon spectrum in QCD, this leads to the quadratical dependence of the 
modification of fragmentation functions or the effective parton energy loss on the nuclear size Ra- 

The LPM interference effect in this study is critical for applying the LQS generalized factorization 
to the fragmentation processes. Because of the LPM interference effect, the gluon radiation with 
small transverse momentum is suppressed. Equivalently, the formation time of the gluon radiation 
has to be smaller than the nuclear size Ra- This requires that the radiated gluon have a minimum 
transverse momentum £ T ~ Q 2 /A 1 / 3 . For such a value of transverse momentum, we can still take the 
leading logarithmic approximation £ T <C Q 2 in the study of fragmentation function in a large nucleus 
(A^> 1). Therefore, the applicability of the our study here is limited to experimental situations in 
which Q 2 must be large enough such that Q 2 ^> Q 2 / 'A 1 / 3 ^> Aq CD . In the meantime, the resultant 
nuclear corrections proportional to a s A x / 3 / '£ T ~ a s A 2 / 3 /Q 2 from double parton scattering are still 
leading for large values of Q 2 . This is why the nuclear corrections to the fragmentation function 
and the induced parton energy loss depend quadratically on the nuclear size Ra- In the processes 
with large final transverse momentum £ T ~ Q 2 as studied by LQS |Ts|| , the leading correction is 
instead proportional to a a A x / 3 /Q 2 . This is because one can neglect the interference contribution if 
1% » Q 2 /A 1 ' 3 . 

We have also considered double-quark scattering processes. Though their contributions to the 
QCD evolution equations can be neglected as compared to quark-gluon scattering, they do have 
a leading-order contribution which mixes quark and gluon fragmentation functions. Since they 
involve quark-antiquark correlations in nuclei, the modification to quark and antiquark fragmentation 
functions will be different. This might give different modification to the spectra for negative and 
positive hadrons as observed in experiments j29] 

There is currently little information on the twist-four parton matrix elements in nuclei, especially 
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the off-diagonal ones. We have only outlined a very crude estimate, assuming factorization of the 
two parton correlations in nuclei. This enables us to estimate the nuclear dependence of the final 
results. Future work on this is necessary in order to have any quantitative study of the problem. 
With that information, one should be able to numerically solve the QCD evolution equation for the 
modified fragmentation functions. 

The method we developed here to study the modification of parton fragmentation functions cannot 
yet be applied directly to other processes, such as high-energy pA and AA collisions. One can, 
however, find important implications from the study in this paper. The A 4 / 3 nuclear size dependence 
of the twist-four parton matrix elements relies on the color confinement within the nucleon radius rjv- 
So the strength of parton correlation as represented by the constant C in Eq. ( f)l|) is proportional to 
rjv- If we replace the nuclear target in DIS by a droplet of quark-gluon plasma, then the integration 
over — y2 in Eq. ([)!]) is no longer restricted to but rather to the effective parton correlation 
length 1/fiD with /xd being the Debye screening mass in the quark-gluon plasma. So the modification 
of the parton fragmentation function and the effective parton energy loss will be sensitive to ■ If 
there is a dramatic change in the parton correlation during the QCD phase transition, one should 
then expect a similar behavior in parton energy loss. 
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APPENDIX 

In this Appendix we calculate the hard part of gluon radiation processes associated with multiple 
quark-gluon scattering in DIS off a nucleus. We use the technique of helicity amplitude for high- 
energy parton scattering, where one can neglect the transverse recoil induced by the scattering. The 
final results will agree with the complete calculation in the limit of soft radiation. Such an exercise 
can help us to cross-check the results of our complete calculation. More importantly, it can help us 
to understand the physical processes more clearly on the level of amplitudes and to reorganize the 
final results in this paper according to the nature of different processes. 

A. Single Scattering 

Assuming the dominant component of a fast quark is i q » [0,£l",0j_]. One has, 

u(£ q , X)^u(l' g , A') « 5\\i rf , (100) 

where n = [0,1, Qj_], A and A' are quark helicities. According to the collinear factorization of 
the parton matrix elements, each initial quark line contributes u(p) while each initial gluon line 
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(with Lorentz index a) contributes p a to the amplitude. Each initial parton line with momentum 
Pi — [xiP + , 0, kr] also has a momentum integral dxi/2n with a phase factor cxjp(—ixip + y~ +iltT-yTi)- 
We summarize these special Feynman rules for initial parton lines as: 

dxi 



initial quark — ► u{p) 
initial gluon — > p a 



2tt 



dxi 
~2^ 



-ixip y ■ -\-ikT-yTi 



(101) 



All the internal and final external parton lines follow the normal Feynman rules. 

We will work with an axial gauge A~ = 0. The final gluon with momentum £ is assumed to carry 
1 — z momentum fraction of the struck quark and has polarization e(£), 



1 = 
e{t) = 



2(1 -z)q- 

e*T • It 
(1 - z)q~ 



-,{l-z)q ,£ T 



,0,er 



(102) 



According to the above Feynman rules, the amplitude for gluon radiation in the single scattering 
case of DIS as shown in Fig. |lj(a) is 



dx 



M*(y)= I —u{xp + q) 1 ^u{p)M" 



tS, , „ er-^Tm -i(x B +x L )p+y- 



M (y) = 2g- 



-T,e~ 



where we have used the on-shell condition for the final quark 

{xp + q-l) 2 = 2p+q-{x~x B - X L ) = , 



(103) 



(104) 



and assumed soft radiation limit 1 — z — > 0. T c is the color matrix in the adjunct representation 
with the color index c for the radiated gluon. 

To calculate the hard part of the parton scattering, one can simply take the square of the amplitude 
Mu^\o)M^ (y), average over the initial quark's spin, sum over the final gluon's polarization and 
integrate over the final partons' momenta. One should also average over initial and sum over the final 
partons' color indices. Note that one of integrations over x is carried out by the overall momentum 
conservation 27r£(^a;i). In the other integral dx/2n, we absorb the factor 1/2tt into the definition 
of the parton matrix elements associated with the processes. Following this convention, we have the 
hard part of gluon radiation induced by single scattering, 

6(x - x B ) 



7T (1) 



dx-Tv[p ■ 77 AI (xp + q) ■ ^]2tt- ^ 



d 2 £ T dz 
2(2tt) 3 1 - z 



M S (0)M S \y) 



dxHl°J(x,p,q) 



Ota 

2^ 



d£l 



-dzCp- 



J(x+x L )p+y 



(105) 



One can see that the above is the same as the hard part in Eq. (pl|) in the soft radiation limit z — ► 1. 
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(a) 



t 



xip + k T 

y vi 

(b) 




^\xip + k T 

v vi v m 

(c) (d) 
FIG. 14. Gluon radiation from a single scattering (a) and double scattering (b-d) in DIS. 



B. Double Scattering 

For double scattering processes as shown in Fig. |l4|(b)-(d), there is one loop-integration which is 
carried out by the contour integration around a pole in one of the two propagators. We also assume 
that the initial gluon carries transverse momentum kx- The on-shell requirement of the final quark 
is 



[q + (x + xi)p + k T — C] 2 = 2zp + q (x + Xi — x B - x L — x D ) = , 



(106) 



where xl and xd are defined in Eq. (|36|). The kinematics in the process of Fig. [l4](b) only allows 
one choice of pole. The contour integration around this pole gives 



dx e~ ixp+ ^ y ~ Vl )^^ x b+ x l+x D )p + v t 



2tt 



{xp + q) 2 + ie 



2p+q- 



-ix B p + y -i{x L +x D )p + y 1 Q(y- —y~\ 



(107) 



The momentum fraction carried by the initial gluon is then x\ = Xl + xd. Since this momentum 
fraction is finite when — 0, the secondary scattering is considered hard. The final gluon in this 
process is from the final state radiation of the secondary scattering. This is what we refer to as 



double hard scattering. Using the helicity selection rule in Eq. (100), one can find the amplitude for 
the process in Fig. 0(b), 



M D(1 \y, yi ) = 2g^ T ^T c T ai e- l< - XB+x ^ p+ y~- lXDP+ yi 
x e ix ^ + (v--y;)i g e(y- -„-), 



(108) 



where a\ is the color index of the initial gluon. Note that there should also be an integral over kr- 
We have absorbed this together with the transverse phase factor exp(ifcr • yrr) into the definition 
of parton matrix elements. 

In the other two diagrams of double scattering in Fig. |lj, there are two possible choices of poles. 
One choice of poles corresponds to soft rescattering of the final quark [Fig. |lj(c)] or of the final 
gluon [Fig. |lj(d)] where the initial gluon's momentum fraction x\ = xd [for Fig. H(c)] or x\ = 



35 



—zxd /(I — z) [for Fig. |l4|(d)] goes to zero for kx — 0. The final gluon is induced by the hard photon- 
quark scattering. This is what we call hard-soft double scattering. The second choice however gives 
the initial gluon finite momentum x\ = xl + xd and the gluon radiation is from the initial state 
radiation of the hard secondary scattering in what we call double-hard scattering. One can find the 
amplitudes of these two diagrams as 

M D{2 \y, yi ) = 2g^ Z T ai T c e-^ XB+x ^ p+y ~- tXDP+ ^ 

x [1 _ e"^ +( -y--y^}ige(y- - yD , (109) 

M D{3 \y,yi) = -2g eT 'S tT Z kT \ T ail T c ]e-^ XB+x ^ p+y ~ ~ lXDP+y i 
(£ T - k T ) 2 

x { e -ixn P + (y--y^)/(i-z) _ e i XLP +( y -- y -)^ ge ^ y - _ y -j _ ^ 110 j 

Since the gluon radiation induced by the hard rescattering in these two diagrams is initial state 
radiation, the amplitude has the opposite sign with respect to the hard-soft processes where the 
gluon is produced by final state radiation of the first hard scattering. 

C. Triple Scattering 

For gluon radiation associated with triple scattering, there are all together 7 different diagrams 
as shown in Fig. [l5]. we list the amplitudes here: 

M T(1) (y,yi,ya) = 2g £ -^^T c T a2 T ai e-^ XB + x ^ + ^ '-^vHvl "-ife) 
M T(2) (y,2/i,y 2 ) = 2g e -^^T a2 T c T ai e- l( > XB+x ^ p+ y~-* XDP+ ^-y^ 

x (-<? 2 )%r -yMy-- vi) , (H2) 



M n3 \y,y u y 2 ) = 2g e -^^T a2 T ai T c e^ XB+x ^ + y~- lx °v + ^- y ^ 

x [1 - ^ x ^ + ^-y^]{-g 2 )e{y- - y-)9(y- - y-) , (113) 



M nA \y,y uV2 ) = 2 3 ^^M[ TQ2;Tc ] Taie -»(^+^) P + y--«^ + fer-^) 

(£ T - M 

x [e-'^D -xd)p + (vi -y 2 )+ixLP + (y~-y 2 ) 



e 



(t- z /(l-z))x DP +( yi -y 2 )-Hxz,P + (y - Vl )](-g*)$(y- - y-)0(y- - y-) , (114) 



M T(5 W,y 2 ) = 2g " T ' {iT =, k ^ T a2 [T ai ,T c }e-< XB+x ^ + y--^ + ^-y^ 



it — kr) 2 



x [ e «z.P + (v--vT) - e - iXDP+ ^~-^y^](-g 2 )e(y^ - y 2 )9{y- - yf) , (115) 



3G 



M Ti6 \y, yi ,y 2 ) = 2g 



€t {It kx) r ,_ i ( XB+XL ) p +y -ix D p + {y -y ) 

(£ T - k T y 



D ix L p + (y -y x ) _ -ix D p + (y -y 1 )/(l-z) 



x (-g*)e(y- -y-)6(y- -Vi). 



J(l-z/(l-z))x D p + (y 1 -y 2 ) 



(116) 



M T{7) (y, yi ,y 2 ) = 2g £ -^^{T ai ,lT a2 ,T c }}e-^ B+x ^y~- lx °P + ^-y^ 
x [1 — e iXLP+ ^ v ~~ v ^^]e i ^ 1 ~ z ^ 1 ~ z ^ XDP+ ^~ v ^^ 

x (-. 9 2 )%r 2/r). 



(117) 



In M , M and M , where the radiated gluon is coupled to the second initial gluon through 
a triple-gluon vertex, we have made variable change Ut — > — fcy. 



xp 



2/ 2/i 



2/2 



Xp 



ojxiP + k T o|x 2 p + k T 
V 2/i 2/2 (b) 



xp 



2/ 2/i 



2/2 



xp 



ZlP + ijoj x 2 p + k T f. 
I 2/1 t/ 2 (d) 




Xp 



xp 



' Xip + fcj 

2/ 

7 jJ 


o x 2 p + £; T of 

2/1 2/2 


. xip + k^ x 2 p + kj 
2/ 2/i 


2/2 


2 q> 


*° £ £ 


' xip + kj 


o x 2 p + &ifo 



xp 

^ 2/i 2/2 (g) 

FIG. 15. Gluon radiation from triple scattering in DIS. 
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To obtain these amplitudes, two loop-integrations are carried out by contour integrations around 
two poles. For each diagram, there are sometimes two possible choices of the two poles, and con- 
sequently two contributions to the amplitude. They represent soft or hard rescattering processes. 
However, one of the quark-gluon scatterings must be soft. In Fig. |l5|(a), for example, the gluon 
induced by the third (hard) scattering after the quark has a soft second scattering. Fig |l5|(b) has 
two contributions. One corresponds to gluon radiation induced by the third scattering, the other 
by the second scattering. One of the two quark-gluon scatterings must be soft. In Fig |l5|(c), one 
contribution comes from gluon radiation induced by the second scattering followed by soft quark- 
gluon scattering. The other contribution corresponds to gluon radiation from the first photon-quark 
hard scattering and double soft scattering afterwards. Such analyses of each diagram help us to 
reorganize the amplitude according to the physical processes. 



D. Soft and Hard Rescattering 



We can reorganize all the radiation amplitudes of double and triple scattering according to our 
classification of soft and hard rescattering. The total amplitude for the hard-soft double scattering 
in Figs. |ii|(b)-(d), as represented by Fig. [?], is 



x 2g 



ct-It e T -(l T -fcx) r r rl ixD p+{y-- v -)/{l-z) 

y2 - L ai- L c - r ., 



{It ~ k T ) 2 

The amplitude for double hard scattering, as represented by Fig. pi is 



(118) 



M%, yi ) = e -^B+^)p + y~-™°p + yi ig e(y- - y -) 



M T cT ai ] 



ex - It £*t • (^z 



{It - k T f 



(119) 



In triple scattering processes, the gluon radiation can be induced by the initial photon-quark hard 
scattering which is then followed by two soft scatterings. We denote the gluon radiation together 
with the first soft quark-gluon or gluon-gluon scattering by the effective diagram in Fig. 0. Similarly, 
the second soft scattering can be either quark-gluon or gluon-gluon as shown in Fig. |l(](d) and (e). 
The amplitudes for these two diagrams with double soft scattering are, 



M T s{1) {y, yi ,y 2 ) = e ^B+^)p+y--™°p + iy;-y-. ) { - g i)o { y- y -)e( y - - y ~) 

x 2g 



Z T-lT £T-{^T-k T ) r T] tXD p+(y--y-)/(l- z ) 



{It - k T ) 2 



(120) 



M T s{2) (y,y 1: y 2 ) = e ~^+^ + ^ -«"P + («r -»,") ( -^ )%r _ y - )e{y - _ y - } 
e T ■ (£ T - k T ) r 



x 2g 



(h 



k T ) 2 



-T a \T a T c ]e~ iXDP ( y ~ Vl 



+ 



er • £t 



[Tax i P"eJ2 j T c 



J(l-z/(l-z))x D p+{ Vl -y 2 ) 



(121) 



38 



where M s ^ is the combination of the first term in M and the second term in M ' while 
M 5(2) i s the combination of the second term of M and the first term of M 

In hard rescattering, the gluon radiation is induced either by the second or the third scattering. 
The amplitude for induced radiation by the third scattering as shown in Fig. |l^(c) is 

M T H{1) (y, yi ,y 2 ) = e-^z+^ + v--^ + (vi-y;){~)g*e{y- -y-)9{y- -y^) 



2g[T Cl T a2 ]T ai 



(It - k T f 



CLP + (y -y 2 ) 



(122) 



It is the combination of M , the second term in M , and the first term in M . If the gluon is 
induced by the second scattering, the third soft scattering can either be quark-gluon or gluon-gluon 
as shown in Figs. |lO](a) and (b). Their amplitudes are 



e*T - It e*T • (£t — &t) 



(It - k T f 



Jx L p+(y -y 1 ) 



(123) 



M 



H(3) 



(y,y 1 ,y 2 ) = e -^+^^ + ^-^ + ^-^)(-). 9 2 %r-y 2 -)0( 2 /-- 2 /r) 



2g[T ai ,[T c ,T a2 ]} 



€t - It e*T ■ C&r ~ ) 
4 (It -k T ) 2 



ix L p+(y -y x ) 



x e 



i(l-z/(l-z)):r. D p + (i/ 1 -y 2 ) 



(124) 



respectively. Again, M H / 2 ^ is the combination of the first term in M , the second term in M 

1(3) 



and the first term of M . M H ^ is the sum of the second term in M , the first term of M 



T(7) 

and the second term of M . 

With the amplitudes of the re-classified diagrams, we can easily calculate the partonic hard part 
of double scattering, 



H 



■D 



d 2 l T dz 



^??3 t^[ mD (°' v^ M °\y^ wi) + mT (^ w> yiW s \y) + M s (o)M T \ y , yi ,y 2 )\ , 

2(27r) J 1 — z 

(125) 



which includes both the double scattering and the interferences of triple and single scattering. The 
final results are the same as listed in Sec. [II B with the splitting functions replaced by their form 
in the limit of z — > 1 . 
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